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ABSTRACT 



Long term numerical integrations were performed with a two- level 
model utilizing the quasi-geos trophic equation set. Friction was in- 
corporated through the vertical derivative of the eddy stress and heat- 
ing was applied as a linear function of latitude. The long-term 
interactions between one or two waves in the zonal direction and the 
mean flow were examined. 

The stability of the initial mean flow was investigated with lin- 
earized equations and the most unstable wave numbers were determined. 
The nonlinear equations were integrated using the most unstable wave 
in the initial conditions. 300 day forecasts were made with various 
values of the heating, |3, wall separation, and wave number. A second 
wave number was introduced into the experiments in order to determine 
the most likely wave number for the fully evolved mean flow. Constant 
amplitude, baroclinic, propagating disturbances were obtained in every 
case where the wall separation was 4,000 km. The experiments with the 
wall separation of 8,000 km produced non-steady waves with fluctuating 
baroclinic and barotropic interactions. 
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I. INTRODUCTION 



In a classic experiment, Phillips (1956) performed a long-period 
numerical forecast with a two-level quasi-geostrophic model in a zonal - 
ly periodic channel on a Beta plane. Friction and heating were includ- 
ed in his model with the latter as a linear function of latitude. 

Although Phillips' experiment did not reach statistical equilibrium, 
he demonstrated that the principal features of the general circulation 
of the atmosphere could be reproduced with a simple numerical model. 

More refined studies involving the primitive equations have been car- 
ried out by Smagorinsky (1963), Mintz (1965), Manabe and Smagorinsky 
(1967), Manabe (1969), and Kasahara and Washington (1967) and (1971). 
These studies which more accurately model the actual atmosphere are 
sufficiently complicated so that it is difficult to understand the 
basic interactions between* the mean flow and the disturbances. 

For this reason, we wish to re-examine Phillips' original simple 
model. In this study, Phillips’ model is further simplified to repre- 
sent the x-variation of the disturbance with one or two Fourier compon- 
ents. This restriction is partially justified by the fact that one 
predominate wave number was observed. Similarily, in dish-pan experi- 
ments with inner cores, single wave numbers were observed in many in- 
stances; see Fultz (1959) and Hide (1953). Herring (1963) and (1964) 
has used the technique of examining the interaction between a disturb- 
ance of sinusoidal variation and a mean field in his study of convection 
between two rigid plates. His estimates of the heat fluxes between the 
plates were within 20$ of the value obtained from laboratory experiments. 
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This suggests that a simplified dynamic model of the atmosphere which 
involves the interaction between a simplified disturbance and a mean 
flow can give reasonable approximations of the northward fluxes of heat 
and momentum. 

Therefore, in this study, the complexities have been kept to a 
minimum to simplify interpretation of the interactions. The two-level 
model utilizes the quasi-rgeos trophic equation set. Fluid motion in the 
model is restricted to the mean flow and one or two wave numbers. The 
mean flow and Fourier amplitudes of the two disturbances are allowed to 
vary in time and y-space; hence an adequate description of barotropic 
as well as baroclinic interactions are possible. 

However, such a drastically simplified model must omit many details 
of the real atmosphere. In the case of this model, a number of physical 
properties of the earth's atmosphere are missing. Two waves in juxta- 
position will not interact directly, but can only influence each other 
through the modification of the mean flow. No oceans, moisture, moun- 
tains, or east-west heating gradient are allowed in the model. Conse- 
quently, this model appears to be one of the simplest models which can 
describe some of the principal mechanisms of the general circulation. 
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II. THE MODEL AND FORECAST EQUATIONS 



The model used in these experiments utilizes the simple two-level, 
quasi-geostrophic equation set. The atmosphere is divided into four 
layers of thickness Ap/2, numbered 0 to 4 in Figure 1. 



0 
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<D 



0 
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AP 



60 



2 

3 



2AP 



60=0 
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Fig. 1. Two-level model. 

The earth's surface is assumed to be a smooth, flat, land surface 
with no oceans. Therefore, the vertical motion at the top and "bottom 
of the atmosphere is taken to be zero. 

Internal friction is represented by the vertical derivative of the 
eddy stress T . 

The quasi-geostrophic vorticity and thermodynamic equations are 



+ W • V£ + P q v - fjj = - Ik • V x g 

|4 +W . v |4 +(Dff ..W 

ot op op p 



( 2 . 1 ) 



( 2 . 2 ) 



where fi = R/C p , S = net heating from radiation, and a, the static 

stability, is 
n 2= 



a 55 ^ [y A - 7] = [ C“ + if ] = 2.3mV 2 cb" 2 

P g P 



= R_T 

2 L/ d 
P g 



(2.3) 



T is taken from a standard atmosphere. 
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Apply Equation (2.1) at levels 1 and 3, which yields 



ST 

*3 

ST 



+ \v„ 



+ vl 



o 1 



o 3 



f ^ 
o Ap 


= Jlk -v X (t 2 - t o ) 


(2.4) 




= ^llc • V x(t 4 - t 2 ) 




f -2 
o Ap 


(2.5) 



If we apply Equation (2.2) at level 2, we obtain 

- 






ip>2 + W 2 



V( 5?>2 + “2 a = 



Ap 



( 2 . 6 ) 



The heating function, S^ 9 is restricted to a linear function in y 
as utilized by Phillips (1956). The model has heating in the south and 
cooling in the north. 

S 2 H h(| - y)/w (2.7) 

where w is the north-south wall separation distance. 

The amount of net heating is proportional to h which varies from 
0.001 to 0.016 kilojoules ton ^ sec 

The eddy stress T q at the top of the atmosphere is taken to be zero. 
The eddy stress term T ^ at level 2, takes the form 

2 aw 

p sr 



T 2 S 



( 2 . 8 ) 

where y is the eddy viscosity. If we follow Smagorinsky (1963), the 

finite difference form becomes 

t 2 = - (P7 ) 2 (SP 2 /Ap> ('V 3 - \\) (2.9) 

where (p g/Ap) = 1/7.9 km ^ is the inverse depth of the 750-250mb layer 
50) _ 4 

x 10 ton m s is an exchange coefficient. The 



and (py) 2 = 



225 

500 



smallest estimates of (p y) ^ was proposed by Rossby and Montgomery (1935) 
for stable conditions, while the largest was utilized by Riehl (1951). 
The intermediate value, 225, was proposed by Palmen (1955). 
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Thus, the stress term at level 2 becomes 



where 



T 2 " - C 2 < W 3 ' V 

225 j x 10 6/79 ton s ^ 
500 ) 



( 2 . 10 ) 

( 2 . 11 ) 



The majority of the experiments were conducted using Palmen's exchange 
coefficient of 225. 

The eddy stress term at the surface, level 4, as formulated by 
Phillips (1956) is assumed to be in the direction of the surface wind 
and proportional to the square of the wind speed. 

T 4 S 2 °D KK (2 - 12) 



where is the drag coefficient = 0.015. 

The surface geostrophic vorticity is approximated by the 750mb 
vorticity 



£4" 





thus W. = W . 

4 3 

So, by defining | Cp |\V^ | 

T 4 " C 4 W 3 

The average value of the 750mb wind speed was taken to be 
Thus, our assumptions make a linear function of wind speed 
= 3.69 x 10 ^ ton s ^ 



(2.13) 

(2.14) 

, -1 

5ms 

with 

(2.15) 



Therefore, as boundary conditions at the top, internal, and bottom 
layers of the atmosphere 

CD = CD, = 0 

o 4 

T = 0 
o 

T 2 = * C 2 <W 3 - V 
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( 2 . 16 ) 



T 4 = C 4 W 3 



For convenience, the following quantities are defined 

* , 

2 

* . iii 

Y t 2 



(2.17) 



implies 



*1 = i + <t> T 

^1 * ^3 






^1 + 



^3 ^1 “ 






^1 " ^3 



- C T 



the mean quantities may be interpreted as applying at level 2. 

d(l> 

Integrate the hydrostatic equation = " a from level 1 to 3, 
which yields 



i - f m ( W 



(2.18) 



Equation (2.18) shows that ^ is proportional to the mean temperature 
of the layer. 

If (2.16) and (2.17) are introduced into the thermodynamic equation 
(2.6), we obtain 



5t *r + \ ‘ - 



o^aAp ks 2 



(2.19) 



The mean vorticity equation is obtained by adding Equation (2.4) 
and (2.5) and dividing by two. 

It V \ * VS M + W T • + P o V M = ilr ,k * V X T 4 (2 - 20) 

The thermal vorticity equation is obtained by subtracting Equation 
(2.4) from (2.5) and dividing by two. 
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( 2 . 21 ) 



+ V v Si + w x • % + e o v T - f Q 4 - 

- f^lk • V X T 2 + f^lk • v X T 4 

Solve (2.19) for and apply this solution to the thermal vorti- 
city equation (2.21). 

Substitute the following geostrophic relations 

W = i-lk x V {> , C - J- V 2 b 



into Equation (2.20) and multiply by f Q , which becomes the Mean Vorti- 
city Equation: 

It v2 ^k + £- ( lk x V( ^) • v < v \) + f (Ik X v<i> T ) -v<v 2 4 ^) 
o o 

d<k 

+ - k 4 - v2 i } = 0 (2 - 22) 

Introduce the geostrophic relations for W and Q into (2.21) and 
multiply by f Q , which gives the Thermal Vorticity Equation: 



where 



|^(V 2 - I4 2 ) + i (Ik x V^) • V(V 2 - H 2 ) ^ 

O 

+ ? ( Ik x v<t> T ) • V(V 2 j^) - 5^ + K 4 (V 2 ^ - v\) 

O 



? ='V 2 *r 


(2.23) 


2f o , in -12 -2 

j - 3.48 x 10 m 


(2.24) 


(Ap) a 




C 4 S -6 -1 

zEf = 1 - 12 x 10 s 


(2.25) 


~4~ - 3 * 59 x 10 s 


(2.26) 
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It can be seen that the frictional coefficient at the surface is 



roughly three times as large as the internal coefficient. 

The mean and thermal vorticity equations (2.22) and (2.23) are the 
prediction equations in the two unknowns (j^ and ((^. . At this point, 
the geopotential height fields are restricted as follows: 

^ = E(y,t) + A(y,t) cos kx + B(y,t) sin kx 

+ G(y,t) cos mx + H(y,t) sin mx (2.27) 

tip - F(y,t) + C(y, t) cos kx + D(y,t) sin kx 

+ P(y,t) cos mx + Q(y,t) sin mx (2.28) 

where E(y,t) and F(y,t) are the amplitudes of the zonal mean flow; 
A(y,t), B(y,t), C(y,t), and D(y,t) are the Fourier amplitudes of the 
disturbance of wave number k, and G(y,t), H(y,t), P(y,t), and Q(y,t) 
are the Fourier amplitudes of the disturbance of wave number m. Both 
waves are in the x-direction only. 

Since E and F are functions of y, barotropic instability is possible 
in this system as well as baroclinic instability. If the definitions 
of ^ and (J).^ are substituted into W = -j-[k x V(|> , we obtain: 
i i 



"m = - f ST 

O J 

l 3 e 
"m = * f 57 

i°**r 
"t 7 5T 

o 

1 ^F 
U T " " f o §7 



f c)x 
O 

1 y T 



v = 1 |£-0 

T j dx 

o 



(2.29) 

(2.30) 

(2.31) 

(2.32) 



where the bar ( ) indicates the x-average. 

At the north and south boundaries, the following conditions are 
imposed initially and for all time: 
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at y = 0, y = W 



“m = “t = ”H = “t = 0 
9u m 9u t „ 

sr = sr = ar = sr = 0 



at y = 0, y = W 



where W is the northern boundary. 

Substitute (2.27) and (2.28) into the mean vorticity (2.22) and 
thermal vorticity (2.23) equations. Now if we separate the various 
terms, and neglect all terms of wave number 2k and 2m and combinations 
of m and k, we obtain ten equations in ten unknowns. 

When we equate the coefficients of the cos kx terms of (2.22), we 
obtain: 



a .a 2 . 2, , k f a 2 B se , a 2 D sf s 3 e 

5t (- 2 -k )A = j t-2 ^ + -2 §7 - B-^ 
ay o oy oy ay 

>2 9 

- P o Bk - K 4 (— 2 - k ) (A - C) 



J5 3 F , 2 ,_c)E , 0F N1 

^3 - k (B §7 + V 3 



y 

(2.33) 



The sin kx terms of (2.22) give: 



a ,a 2 ,2.„ k . a z A se a^c aF , ,a J E .^f , . 2 .. as , .ar. , 

■ V'Z? 5? ■ * * "57 + V 



0 2 a £e 



d 2 C 0F 



a S 3 E 



^ 3 f 



Sy 

+ h* - V^2 

0y 



0y 



0y dy 



k ) (B - D) 



(2.34) 



The cos kx terms of (2.23) give: 



s , k 2 2 _ k r aE a^D aF a^ _ 2 

at V M } £ 0 Sy a y 2 Sy a y 2 



a 3 s. 



a 3 F 

D - 2-4b] - 0 Dk + K. ( 



ay ay 

+ K 2 (Ck 2 



a 2 c . 



0y 2 






-k 2 ) (A-C) 



M 2 ) 



(2.35) 



The sin kx terms of (2.23) give: 



aE a 2 c 



a ,b* ,2 2 . _ k , „„ „ o 

5t (- 2 -k •+* )D = 1 1 - ^ — j 
oy o oy 



dF 0 2 A 






0F a 



= 7 t “^““§7~ + + [i ) + ^;A(k - fa ) 



0y 



0y 
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(2.36) 



3 3 e„ 



3 3 f, 



a 3 



5t 



3y 2E " 



P Ck 4 
o 


k r 2 

K ‘ V 


,-k 2 ) (B-D) 


*h 


: of x 


in (2.22) give: 




A 


° 0 


+ C ~2- 




3y 2 


dy 


3y 2 


3y 


A 


- H— — ^ 


+ ^i- 


3fp, 

Q— d 


by 


dy 


by 


3y 2 



$> 



(2.37) 



k B 

2f d} 
o J 

, m c) 

+ F 5j 

O J 

„ 3 ,3e 3f. 

4 5y by §y 

The cos mx terms of (2.22) give: 

d ,b 2 2. m r ^ 2 H bE , b^O dF u d 3 E ^ 3 F 2 /TT dE , -df\ , 

*<‘5? ’ • ,G ■ V^7 5 s - V -^3 - <% * V> 

^2 

- P ftn - K 4 (— 2 - in ) (G - P) (2.38) 

3y 

The sin mx terms of (2.22) give: 



S A 



5t<7-2- n)H "? [ -rl^-^2^ + G : 

gy r\ r^\r * ri^r * t 



b 2 G bE d 2 P dF , ^ 3 E , t 3 3 F , 2 . _^E , J^F 



dy 



^y 



by by 



+ P o Gm - K 4 (~ 2 - m Z ) (H 
dy 

The cos mx terms of (2.23) give: 

^2 0 „ _ ^2 



Q) 



(G §f + % )] 

(2.39) 



3 ,3 2 2. „ m ,3e 3~o , 3f 3 h 3e_. 2, ,,2. 3f . 2 ,,2. 

?t (- 2 - m - H )P = + ^(m + n )- ^H(m - M ) 

3 3 2 2 

- ^-|q - ^H] - P o Qm + K 4 (— 2 - m 2 )(G-P)+ K 2 (Pm 2 - ^) (2.40) 

by by by by 

The sin rax terms of (2.23) give: 

3 .3 2 2 2. m . 3 e 3 2 p 3f 3 2 g ^ Se,,. 2, 2. 3i* 2 2, 

5t (-2-”> - M )Q = f Q t- §7 ^2 - '37 ^2 + 3^ P(m + ^ > + 3^ (ffi - ^ > 

3 3 2 2 

+ — ?G + ^-§P] + 6 Pm + K (— - m 2 )(H-Q)+ K (Qm 2 - $-$) (2.41) 

'sO'nO O 4 ^ Z Z 

oy ay ay oy 
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The coefficients independent of x in (2.23) give: 
\2 



d , ,2. k S . S 2 D S‘C „S‘b . . .2 .... 

5 t ( r 2 ‘ 4 )F = 2T ?y [A 7T - + c-^ - + la (BC - AD)] 



„d 2 c 



„5 2 b 



,S 2 a 



Sy 



by by by by 



+ w - it + A 

2f o Sy Sy 2 S y 2 dy 2 



O^G 

dy 



f + M- (HP - GQ)] 



d ,dE b¥ 



2 * S 2 



+ K 4 dy ( o^ ' 5^ ‘ M 2 



- K, 



b 2 ¥ 
2 ^y 2 



(2.42) 



Equations (2.33) - (2.42) constitute the prediction formulas. 
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III. PROCEDURE AND GOALS 



The north-south grid distance Ay was set at 200km. A time step. 

At, of 1/2 hour was chosen to allow wind speeds up to 100 m s ^ with 
computational stability. 

Centered time differences were used for all quantities except those 
involving friction. The frictional terms were computed at time t - At. 
The first time step was a forward step. The Euler -backward time step 
was introduced every 24 hours to selectively damp high frequency waves. 
These waves are produced by the separation of solutions at even and odd 
time steps caused by the leap-frog time scheme. 

The finite difference scheme used is illustrated below with the 
function F, 



| = 2fe (F j + l * F j-1> 

dy 2 (Ay) 2 <Fj+1 ^ + - 

57 ' 2P J«* V" P J ' !P J-i + P i- !>l 

where j is the grid index and Ay the distance between grid points 
(200km) . 



Equations (2.33) through (2.42) except (2.37) were solved by a 
Gauss pivotal elimination matrix method described by Richtmyer (1957). 
Equation (2.37) was solved by a direct marching process). 

The model was run as a Beta-plane experiment with periodic east- 
west motion between north-south walls of separation W, 4,000 or 8,000 
km. 
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The 8,000 km experiments were likened to the atmosphere. It was 

hoped the mean temperature gradient from the north to south boundary 

at level 2 would be comparable to that of the atmosphere at 500 mb. 

The source of available potential energy, (S*P), should yield a value 

not much less than 60$ of the atmosphere’s total zonal generation of 
-4 -2 -1 

56 x 10 kj m s formalized by Dutton and Johnson (1967). It should 
be noted that this model contains no east-west heating gradient; there- 
fore, conversion of heat into available potential energy (APE) will be 
less than observed. 

Lastly, the winds at 250 mb should exhibit a jet structure and the 
750 mb winds portray the belts of zonal easterlies and westerlies. 
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IV, VERTICAL MOTION EQUATIONS 



Vertical motion at 500 mb is calculated through the thermodynamic 



equation (2,19) 
2 



m 2 " A P a 



( . ,dC k /t Ae . r dD k A dF N , 

Jcos kx[^ - - B^)] + sin kx[^ + -^(C^ - A^)] 



r c)P m , A, , . r ^Q , m /t Ae ^F x1 

+ cos mx[^ - - l^)] + sin «k[^ + ^(P§^ - <fy)] 



+ 9F HS 2 

+ 5t ‘ T 



Apaf 



1 T r^A I -n^C , r,^B a SD-. 

k[ -^ + % + c §7 ' V 



, Ag , J)P , Am Aq, ) 

+ m[ - % + % + % - % ] j 



(4.1) 



all terms of wave number 2k or 2m and combinations of m and k have been 
neglected. 

The vertical motion, cd, can be written as 
cd = cd + cd' 

where the average vertical motion is 



c)F 2 



d 



= 55 S-T + F [k Sy (BC - "» + m 5y (PH - «»1 < 4 * 2 > 



again, the barred quantity means x- averaged. 

The disturbance vertical motion, cd' , is given by 

cd * = cd' (cos kx) 4- cd' (sin kx) + cd' (cos mx) + cd' (sin mx) 
c s c s 



where 
cd^ ( cos kx) 

cd' (sin kx) 
s 

cd^ ( cos mx) 

cd' (sin mx) 
- s 



2 I i r^C k ,_dE ^F n1 ) 

55 I kx[ 5t • T ( % ' V 1 ! 

2 \ . . r dD , k ,_dE a ^F n1 1 

55 i sln ^ + T { % ' V ] ) 

55 ( cos • r ( % - V 1 1 
sz? j sininxC §i + f; (I ^** % )] \ 



(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 
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v. COMPUTATION OF AVERAGE NORTH- SOUTH WIND 



The average north-south flow computations coupled with the average 
vertical motion provides one with a fast method of calculating the cir- 
culation in the meridional plane; i.e. Hadley or Ferret cells. 

If we apply the continuity equation at level 1, we obtain 



^ = - V • W. 

op 1 



dv. 






5y~* dp 

Thus, the finite difference form is 



(v,), = (v_) 



l'l 

where ( v ^)^ = 0 



. £ 3 . 



1'J . Ap 



(m,). + (a) ) 



Xl 



lli+i. 



(5.1) 

(5.2) 

(5.3) 
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VI. ENERGY AND ENERGY TRANSFORMATION EQUATIONS 



The energy equations were obtained from Phillips (1956). 
Mean Potential Energy per unit mass. 




Disturbance Potential Energy per unit mass. 




Mean Kinetic Energy per unit mass. 
Disturbance Kinetic Energy per unit mass, 



( 6 . 1 ) 



( 6 . 2 ) 



(6.3) 



C =— 2 f [(V4)’) 2 4- (V{)') 2 ] dy 

4Wf J 1 J 



(6.4) 



If we substitute the (J> T and <j^ into the above energy equations. 



we obtain 

2 N 

P = -i 



2 2 F 2 

2Nf j=0 J 
o J 



U 2 N 2 2 

P' (k) = 2 (C‘ + D‘) 

4Nf j=0 J J 

o J 



2Nf Q j=0 



N 



P' (m) = -“—s 2 (P 2 + Q 2 ) 

4Nf j=0 J J 



. 2 . . 2 , „ 2 , „ 2 , 2 n 



K ' (k)= lo^ + B J + ^ + ^ ] 
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K ' (m) = i -2 j = 0 tm2(G J + h j + p j + Q j } + + 



+ (i^)! + (i 2 ) 2 ] 

dy j dy j 

where the subscripts k and m stand for wave number, and N is the 
number of grid increments such that 

W = NAy 

Following the procedure of Phillips (1956), the energy flow diagram 
of Figure 2 should make the physical interpretation of the energy trans- 




— *( K ^4) 



The flow of energy is in the direction of the arrows if the asso- 
ciated transformation ( • ) is positive. 

The energy transformation equations are: 

lA n 

(S*P) = 2 (S . • F ) 



2Nf Q 2 j=0 j 



iA N 

(P*P')(k) 2 [(A-D - B*C). 

2Nf j=0 3 

o J 



Oj 1 



(P * ’K' ) (k) = 



-1 



2Ap N 



N 

E [C.cn'.(cos kx) + D.<D'.(sin kx) ] 



j=0 



j cj 



j sj 
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,a 2 B . „s 2 c 



,dF x 



^d 2 C a ^ 2 D ^ 3 2 A ^ 2 B n 



S 2 d 






* Vj ‘t? ‘ ' *■*>•' 
1 



N 

2 



9y^ Sy 2 Sy 2 ^y 2/jJ 



(P<K) = VT ^ < F j * ®J> 



(k*t 2 ) = 



N f 



IN >>2 

j-o i^-fv fj 



Sy' 



2 j 



(K'-T 2 )(k) = 



(k*t 4 ) = 



A 



-^- 2 l ! k V + DV(4 + ^l 

2Nf j=0 J 9y 5y J 



N 

2 



(E-F). • - a A> 



N f j=0 
o J 



(K f • T. ) (k) = 



2Nf 



j=0 



V dy 2 j 

-A [(^f - ^§) - k 2 (A-C) ] 

J c)y 2 dy 2 J J 



-k 2 (B-D)J 

J c)y 2 % 2 J J 

The energy transformations of the m wave disturbance have been 
omitted for brevity. By replacing k. A, B, C, and D with m, G, H, P, 
and Q, the entire transformation set will be complete. 
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VII. INITIAL CONDITIONS 



A natural initial state would be one of no motion and constant tem- 
perature. Then the system could run until a steady state is reached or 
some type of instability occurs. However, to save computer time, the 
initial state was chosen to be near a critical state for development of 
the most unstable wave. 

When the mean flow model was allowed to run with no disturbance, 
the wind field reached a steady state; the Reynolds stresses balanced 
the heat input. This steady state was comprised of a meridional Hadley 
cell superimposed on the zonal flow. This Hadley cell mean flow was 
utilized as the initial condition. 

If we set ^- = 0 and A-D, G-Q = 0 , 
r} r}E r}F 

the F equation yields 0 = K 4 £ y (^ - , 

2 

the E equation yields 0 = -(U^tthCr - y) /W - K . 

2 , 2 dy 2 



The solution for F is: 

2 2 3 

F= Oh (Li _ zl- 



a 

k 2 



6W‘ 



(7.1) 



This solution requires that E = F which implies that <|l = (L • In this 



u 3 = 0 and U]L = 2 u t = 2x1^ 



Pedlosky (1964) derived a simple stability criterion for the mean 
flow. If the potential vorticity gradient (g^) changes sign in either 
the y-direction or the vertical instability is possible. 
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For this model bq/by is given by 
dq 1 q ^ U j , “ v 

5f" = p o-^r + F<"i- V 

la! B ^ str -s 

& ^ - 2 < u i - V 



(7.2) 



(7.3) 



Although u^ = 0 initially, easterlies and westerlies develop with 

time as observed in the earth’s atmosphere. However, because u^^ > 

we shall neglect u^ for this discussion. 

If we neglect u in equation (7.3), we see that 
3 2 
p = y-u 

H o 2 1 

is the critical value of the mean flow. Above that value, instability 

may occur for some wave number. If we solve for u, . . ..a value of 

J 1 critical 

^critical ma ^ ^ er ^ ve ^ through the Hadley cell initial condition (7.1). 



In this paper 



The 
A = 200 
B, C, D 
G = 200 
H, P, Q 



h . . , = 0.0011 

critical 



disturbance initial conditions are 
sin TTy/w 
= 0 



sin TTy/w after 100 days 
= 0 
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VIII. RESULTS 



A. HIERARCHY OF UNSTABLE WAVE NUMBERS 

The objective of the first experiment was to find a hierarchy of 
wave numbers in order of wave growth for a particular mean flow. The 
mean flow was given by the Hadley cell solution with h = 0.004 kj ton ^ 
s ^ and it was held constant with time. Since the mean flow did not 
vary with time, the equations for the disturbance became linear. Each 
case was run until expotential growth was achieved. Figure 3 shows the 
graphs of wave number vs wave growth rate for two sets of conditions: 

{3 =0 and |3 = 1.67 x 10 and W = 4,000 km and W = 8,000 km. The 

growth rate shown is 

- . Amp(t ) ? 

growth rate = In [ Amp(t; - ) ] x 10 

where Amp(t) = amplitude of the disturbance at time t. 

Pedlosky (1964) has computed the growth rate for this wind profile 
using different values for various parameters. 

Figure 3 exhibits how the Beta term stabilizes the longer waves. 

B. INTERACTIONS BETWEEN WAVE AND MEAN FLOW 

The second set of experiments were devised to test the growth of 
the most unstable amplifying wave, allowing interactions with the mean 

flow, as a function of heating. All experiments were conducted using 

-4 -1 -1 

Palmen's exchange coefficient of 225 x 10 ton m s and 

heat = h . . . (0.001 kj ton ^ s *) as the initial condition for the 

critical 

Hadley cell. 
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Each experiment was run for 300 days. The k wave was introduced 
immediately; the m wave (if existed) was introduced at 100 days. 

The following tables were set up to briefly describe the character- 
istics of each wave number by way of symbolic nomenclature. The sym- 
bols are: 

D = wave dies immediately 
N = a non- steady wave, disturbance energy 

oscillates up and down as the wave interacts 
with the mean flow. 

SW = a steady wave, total disturbance energy 
2 -2 

is = ( ) m s . 

Peak - the peak amount of disturbance 
2 -2 

energy = ( ) m s at day. 

E'(k) vs E 1 (m) - comparison of k and m wave 

disturbance energy over the 200 day period 
when both waves are introduced. 

E'(m) > E'(k) - if che m wave disturbance energy 

did exceed that of the k wave, then the m wave 

2 -2 

disturbance energy, E'(m) = ( ) m s 

on day. 

- 2-1 

S*P = ( ) kj m s is the average amount of conversion of heat 

to APE usually computed after the mean flow has come to some 
degree of equilibrium with the disturbance. If the system has 
not run long enough to reach the equilibrium, then the value is 
given which best represents the average S*P thus far incurred. 
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C. u 1 AND Uy VERSUS HEATING FOR THE STEADY WAVE EXPERIMENTS 

The table below shows the dependence of the zonal flow upon the 
net heating rate for the Steady Wave experiments. 



heat (kj ton ^ s ^ 


exp# 


U T ( m s" 1 ) 


Ui<m 


0.004 


3 


6.9 


14.1 


0.008 


4 


7.1 


15.1 


0.012 


5 


7.2 


15.6 


0.016 


6 


7.2 


16.0 



Table 5.11 



Table 5.11 lists the steady state values of the zonal flow. The 

rather small increases in u T with heat is best explained by the rapid 
JL 

baroclinic response to the wind sheer decaying the unstable zonal flow. 



D. A BAROCLINIC, STEADY WAVE 



Experiment #5 



= 1.67 x 10 
= 4,000 km 



-11 



k =7 
heat = 0.012 



exhibits a 



very interesting feature in this model: a steady, baroclinic wave. 
Figure 4 is a graph of mean kinetic energy, K, mean potential energy, 
P, mean energy, E, and total energy, E, for experiment #5. Initially, 
the mean flow grows rapidly, as shown on Figure 4, with negligible 
disturbance energy. At the 40th day, the zonal mean flow has reached 
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its peak and decreases rapidly while the disturbance begins to grow. 

By the 115th day, the entire atmosphere has reached a steady state. 

Initially, the maximum u^ = 4.9 m s ^ and u^ = 2.4 m s \ At the 
20th day, the maximum u^ = 14.2 m s ^ (henceforth, whenever zonal mean 
flow is mentioned, it will mean the maximum speed in the north-south 
column) and u T = 7.5 m s \ The 40th day exhibited the maximum mean 
flow energy, due to the maximum P, with u^ = 21.7 m s ^ and u T = 11.1 
m s However, the winds continued to increase until the 46th day, 
as attested to by the K graph, with u^ = 25.3 m s ^ and u^ = 10,9 m s ^ 

Figure 5 (E*(k)) portrays how the wave grew beginning at the 32nd 
day. By the 38th day, the wave responds to the strong baroclinic in- 
stability, as shown by the large u T , and grows expotentially. Baro- 
tropic damping (K ' *K) occurs exclusive of any barotropic growth the 
entire cycle of the wave, and exhibits large oscillations during the 
rapid growth of the wave. 

The wave continued to grow by baroclinic instability until the 56th 

day when the wave reached its peak energy with a disturbance energy, E' 
2 -2 

= 123 m s . By this time, the thermal wind had decreased from a peak 
of 11.2 m s ^ on the 42nd day to 7.5 m s * on the 56th day, destroying 
the source of the baroclinic growth. 

On the 38th day, when the wave began its large baroclinic growth, 
barotropic damping begins to reach its greatest value as seen in Fig. 5 
Initially, as the barotropic damping oscillates about its mean of 
1.4x 10~^kj ton^s^,the baroclinic instability is by far dominant. How- 
ever, by the 58th day, barotropic damping plus the loss of eddy kinetic 
energy, K* , by frictional dissipation becomes larger than P ,# K'; hence, 
baroclinic effects must be diminishing. 
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As the barotropic damping reaches steady state, it appears to come 
into equilibrium with baroclinic growth, hence, the steady state of the 
disturbance energy. 

Figure 6 is a graph of the conversion of net heating to available 

potential energy (S - P). The graph shows the rapid rise of S«P to a 

-4 -2 -1 -4 -1 

peak of 16 x 10 kj m s in a vertical column (1.6 x 10 kj ton 

s at 40 days, then its decrease to the steady state value of 8.8 

-4 -2 -1 

x 10 kj m s by the 104th day. 

At the end of the 300th day run, the entire atmosphere has been at 
a steady state equilibrium since the 115th day. The energy conversion 
terms feeding P, P', K, and K' were constant with respect to time, as 
were the energy terms; therefore, the conversion terms should balance 
for each energy term. Figure 7 is a pictorial display of the 300th day 
energy flow after Phillips (1956), the only conversion not balancing 
being the (K*T^). 

In summary, the steady wave is achieved by the conversion of P to 
K' which is ultimately balanced by frictional losses and barotropic 
damping through the action of the Reynolds stresses. 

E. A NON-STEADY WAVE WITH BAROTROPIC EFFECTS 



quite interesting feature, a wave initially amplified by baroclinicity 
which then grew by barotropic effects. 

Figure 8 is a graph of the mean flow energies and total energy for 
experiment #12. Notice the large and small oscillations in each energy 
graph and how the total energy graph is smoother than the mean flow 




iheat = 0.016 



(3 = 1.67 x 10 
w° = 8,000 km 
k =6 




exhibits another 
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energy graph. This is due to the disturbance drawing energy from the 
mean flow and vice versa. So when the disturbance decreases, the mean 
flow increases and vice versa which tends to smooth out the total 
energy graph. 

On the 20th day, the wave grew rapidly by baroclinic effects. At 

that time, u^ = 9.3 m s ^ and u^ = 18.1 m s ^ while on the 30th day, 

u^ = 10.9 m s ^ and u^ = 30.0 m s Since only baro tropic damping was 

occurring and the mean flow increased since the initial 10 days of 

growth, baroclinic effects amplified the wave to its initial peak of 
2 -2 

69.4 m s at 30 days. 

Figure 9 is a graph of the disturbance energy and the baro tropic 
instability. It may be seen that the barotropic growth term feeds or 
damps the disturbance. At day 30, the wave has peaked due to baro- 
clinic growth. From the 31st to the 33rd day, the baroclinic instability 
subsided and the barotropic damping destroyed over half of the wave 

energy and increased the mean flow to u^ = 15.3 m s ^ and u^= 38.2 m 

-1 

s 

The 34th is the first day that barotropic instability occurred and 

2 

although it only lasted for two days, the disturbance peaked at 77 m 

-2 

s . On the 37th day, barotropic damping again diminished the disturb- 
ance until the 40th day. On the 37th day, u^ = 14.1 m s ^ and u^ = 

37.9 m s ^ while on the 40th day, u^, = 17.2 m s ^ and u^= 42.6 m s \ 

To further supplement this discussion, the K graph is given on the 
bottom of Figure 8, without its energy scale. The graph helps to por- 
tray the barotropic damping effect, gain of mean flow energy to the 
disturbance and vice versa. 
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On the 40th day, the barotropic instability causes explosive ampli- 
fication of the wave. However, in the 42nd and 43rd days, barotropic 
damping occurs and explosive baroclinic growth of approximately five 
times the average P 1 *K’ takes over. 

This baroclinic growth can be seen on the small growth of the K 
graph from the initial loss due to barotropic growth from the 40th to 
the 42nd day to the small gain from the 42nd to the 44th day. 

From the 44th to the 53rd day, growth is a combination of baro- 
tropic and baroclinic effects. Possibly the barotropic term is domin- 
ant, as with each barotropic damping, the disturbance graph changes 
slope or declines, hence a greater dependence on barotropic growth. 

The mean flow has declined greatly due to both a baroclinic and baro- 
tropic drain on it, as shown by the 44th day u^ = 14.2 m s ^ and 
u^ = 40.5 m s 1 while on the 53rd day, u T = 11.9 m s *and u^= 23.9 m s 

From the 53rd to the 82nd day, the wave damped for the most part 
due to barotropic effects, although a small amount of baroclinic damp- 
ing does occur. 

The disturbance peaks on the 56th day, due to barotropic and baro- 
clinic growth, with the K graph reaching its minimum on the 55th day. 
Again, baroclinic effects do create growth in the wave as seen when 
barotropic damping is occurring and the disturbance pauses from its 
decay to grow for a day, then continue its decay. 

Figure 10 is included to observe the variation of S*P term in a 
non-steady wave. Figure 11 is included? to show the many oscillations 
of the P*K term (thin line). An interesting phenomenon noted on Figure 
11 is the comparison of the P*K conversion with the K term (heavy line, 
without energy scale) . The two terms tend to counteract each other 
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Fig. 10. Conversion of Net Heating to Available Potential Energy for Exp//12 
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rather than reinforce each other. In the same fashion, K ,# K also 
counteracts the P*K conversion with K'*K on the average dominating 
P*K. Hence with barotropic damping, (K 1 *K > 0) the P*K conversion 
usually becomes negative. The reverse also occurs. 



graph begins when the disturbance and mean flow reach a state of 
equilibrium after the,, initial growth period of the disturbance on the 
mean flow. The graph runs from 80 to 110 days. The easterlies present 
at 250 mb are attributed to the fact that no horizontal diffusion is 
present in the vorticity equation. 

F. A TWO WAVE EXPERIMENT 



few experiments conducted where E 1 (m) > E'(k). The explanation for the 
dominance of the m wave over the k wave, in opposition to the results 
of section V. A., is attributable to the modification of the mean zonal 
flow by the k wave. Figure 13 is a graph of E'(k) and E' (m) versus 
time. Figure 14 is a graph of the K' *K conversions for each wave. It 
should be noted that, on the average, both conversions act in the same 
direction in quite varying degrees. 

The Pedlosky criteria (7.2 and 7.3) proved to be a necessary condi- 
tion for baroclinic and barotropic instability. The baroclinic criteria 
was met for the initial condition and throughout the forecast. The 
barotropic criteria was met continuously after the 34th day. At no time 
did barotropic growth occur when the Pedlosky criteria was not met. 



Figure 12 is a graph of the u^ flow with respect to time. The 



Experiment #12 




P = 1.67 x 10' 11 
w° = 8,000 km 
k =6 
m =7 

y heat = 0.012 




9 



is one of the 
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Between the 162nd and 164th day, the E'(m) becomes larger than the 
E'(k) and remains such throughout the run. At this time, both waves 
are experiencing barotropic instability. However, the m wave is grow- 
ing and the k wave is decaying. 

The time frame, from 162 to 164 days, is explored further on 
Figures 15, 16 and 17. Figure 15 and 16 is a graph of P' *K' (bold line) 
and K’ *K(thin line) for waves k and m respectively. A unique differ- 
ence is apparent between the two waves. The baroclinic conversion 
P'^K’ never approaches zero with the m wave but does many times with 
the k wave, and even becomes negative. When the baroclinic term be- 
comes negative, baroclinic damping draws from K'. Usually, the P' *K' 
conversion was positive, thus the waves were experiencing baroclinic 
growth. On day 164, both waves are experiencing barotropic growth; 
however, the k wave is undergoing baroclinic decay and the m wave is 
experiencing baroclinic growth. The energy flow chart is shown below 
to help explain the growth of the m wave and decay of the k wave. 



m 

i 

1.41 

J 

| P' |< — 0-16 — | K' 1— *0.38 

l 

026 
Ic wave 




0-44 



032 



m wave 



net change to K'(k) net change to K f (m) 

= +.61kj ton ^ s ^ = +.92kj ton ^ s ^ 

Possibly this difference in net inflow to the K' term for each wave 
gives the growth advantage to the m wave for a ten day period. 
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If we look at the waves on the 164th day, we find an interesting 
feature of the baro tropic growth. Although the k wave K' *K is largest 
(negative), the m wave is experiencing the growth through baroclinic 
instability. Since the conversion terms are the meridional average, 
it is possible the barotropic growth of the m wave in a sector is sig- 
nificantly greater than barotropic growth of the k wave in a similar 
region. This possibility is explored in Figure 17, a graph of the 
phase angle of the k and m waves and u^. The wave is experiencing 
barotropic growth if the tilt of the wave is opposite the shear of the 
mean flow. Figure 17 depicts the k wave growing barotropically south 
of the zonal maximum of u^ to the easterly jet, and the m wave growing 
barotropically north of the zonal maximum to about y = 6,000 km, thus 
demonstrating the greater barotropic instability of the k wave. 
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Disturbance Energies of k wave and m wave for Exp //26. 
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wave & m wave vs Time for Exp//26. 
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Fig. 15. Baroclinic and Barotropic Conversions of k wave vs Time for Exp 1126 . 
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Thousands of kilometers 




Fig. 17. 



Phase Angle and U vs y for Exp #26. 
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IX. CONCLUSIONS 



In order to facilitate physical interpretation of the results, the 
model was restricted to a mean flow and one or two waves. A steady 
state solution which was independent of x was obtained. This flow was 
of the Hadley type with zero mean velocity at level 3 and a parabolic 
wind profile at level 1. The stability of this mean flow was investi- 
gated by performing initial value integrations with the linearized 
frictionless equations. The calculations were continued until expon- 
ential wave growth was achieved. At this point, the growth rates for 
a range of wave numbers were calculated. The growth rate of the waves 
was determined as a function of wall separation and P. The growth rate 
was reduced by the P effect and increased greatly by large wall separa- 
tions . 

The nonlinear frictional equations were integrated from an initial 
state which consisted of the Hadley mean flow upon which was super- 
imposed the wave number of maximum growth rate. For each combination of 
heating, P, and wall separation, a 300 day prediction was made with the 
optimum wave number from the linear study previously mentioned. Each 
of these forecasts was repeated but at 100 days a new wave was intro- 
duced. The purpose here was to determine the largest amplitude wave 
which would exist on the mean flow, the latter having been modified by 
nonlinear interactions by the first wave. A total of 45 experiments 
were made as a function of wall separation, heating, P, and wave number. 

Constant amplitude propagating disturbances were obtained in every 
case where the wall separation was 4,000 km. These experiments 
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contained cases where j3 = 0 and P ^ 0 as well as a large range in heat- 
ing rates. The steady waves gained energy through baroclinic conver- 
sion processes and lost energy through barotropic decay and friction. 

The remaining experiments with the wall separation of 8,000 km pro- 
duced non-steady waves with fluctuating baroclinic and barotropic in- 
teractions. The roles of the different types of instability in the 
growth and decay of the waves were easily distinguished. When one in- 
stability was causing strong growth, the other instability was usually 
a minimum growth conversion or was decaying the disturbance. The baro- 
tropic and baroclinic conversions tended to be out of phase with one 
another, one a maximum - the other a minimum. The P*K conversion, 
however, tended to be a constraint on the K term and adjusted to the 
barotropic term. With strong barotropic damping, the P # K was negative, 
while with strong barotropic growth, the P*K conversion was positive. 
Hence with strong barotropic decay, the Ferrel cell became dominant over 
the Hadley cell. 

Over the long term mean, the waves obtained a net gain of energy 
through P'*K' and a net loss of energy through K' *K. The mean tempera- 
ture gradient from the south wall to the north wall and the S*P conver- 
sion were directly related to the heating and wall separation. With 
the largest heating and wall separation, the temperature gradients were 
somewhat greater than the north-south 500 mb temperature gradients in 
the earth's atmosphere. This is related to the larger than normal 
heating used to offset the lack of east-west heating gradient. Even 
with the largest heating and wall separation, the conversion of heat to 

available potential energy, S*P, is only about 65$ of the average for 
-4 -2 -1 

our atmosphere of 56 x 10 kj m s as estimated by Dutton and 
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Johnson (1967). This underestimate is also the result of the lack of 
east-west heating. 

The zonal mean flow at level 1 exhibited a definite westerly jet 
in the mid- latitudes, which sometimes split into a double jet. Due to 
the lack of a horizontal diffusion term, the horizontal gradients were 
sometimes quite large and occasionally an easterly jet would appear for 
a short time. The zonal mean flow at level 3 exhibited the belts of 
zonal easterlies and westerlies typical of the earth’s atmosphere. 

In this study, we found with a wall separation of 4,000 km only 
steady baroclinic waves, and with a wall separation of 8,000 km 
fluctuating baroclinic and barotropic interactions yield non-steady 
waves. It was determined that the predominate wave number in the fully 
developed flow was in every case close to the value determined by the 
linear stability analysis. In the actual atmosphere, where walls are 
not present, barotropic instability may be important in giving day to 
day fluctuations in the mean flow. 
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